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c3 ■ , Her ^ ' 1 (IGiacconi et al.||1971|) and Cen X-3 (|Tananbaum et al.| as sociated with the phy sical accretion scenario first suggested 

|1972|). The basic mechanisms of the pulsed emission were un- bv lDavidson & Ostrikeil d 19731) . A summary of the model prop- 

derstood quickly: radiation originates from the accretion of ion- erties will be given in Sect.|2] 

ized gas onto a rotating neutron star (NS) from a nearby com- A unique spectral characteristic of many HMXBs is the pres- 

pamon, an O or B star (M > 5M ) in high mass X-ray bina- ence of cyclotron resonant scattering features (CRSFs), which 

nes (HMXBs), a later than type A star (M < M ) in low mass provide a tool for direct measurem ent of the magnetic-field 

X-ray binaries. For strong magnetic fields (B ~ 10 " G) of strengths of accre ting pulsars, a key ingredient of the model. The 

the compact object, the plasma is threaded at several hundred fundamental line appears at E = n 6Bn x (1 + z) -i keV , 

stellar radii and then funneled along the magnetic field lines , n . ,. c \, . ■ c lrv i2,^ 

b b where Btj is the magnetic held strength in units of 10 G 

onto the magnetic poles, forming one or two accretion columns , • .. .. , ■ u e . ■ .? v e 

/o ■ i 9 d imo r\ — ^ — b V> > i Jlim-^ K ti, , and z is th e gravitational redshift in the line-forming region 

d Pnngle& Rees| U97l | Davidson & Ostrike jGpl). The system (Was serma ; & P Shap l^| rj 983) . However, the determination of 

is powered by the conversion of gravitational potential energy 7 I ^ e e A » i r ■ u. u u 

. , . 1 , • , ■ , • r the centroid energy of the fundamental line might be ham- 

mto kinetic energy, which is eventually emitted in the form of ju v i t. j»u» • / 

, & -" , ., , , , ,• • pered by its complex shape, due to photon spawning (see e.g. 

X-rays as the plasma decelerates, possibly through a radiative rt? — — J -r—r . r , i u nn gK t, , , • • f ■ 

i , \ , ^ -r. i • ■ Krevkenbohm et al. 2005). The second harmonic is deeper and 

shock, and settles onto the stellar surface. Pulsation is gener- T ^ rrr , , ^ , £ ^ , 

, . . , ,. . , b , can be described by a simple analytical function such as a 

ated because of the nonalignment between the magnetic and ro- . . .. . f . , .. i XT . , ■ 1 

tational axes Gaussian curve, since it is due to pure absorption ONishimural 

2003). Phase resolved spectral analysis is essential for this kind 
of study because the spectral characteristics of the lines de- 
Send offprint requests to: C. Ferrigno pend heavily on the angle between the observation direction and 
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Abstract 

Context. Highly magnetized pulsars accreting matter in a binary system are bright sources in the X-ray band (0. 1-100 keV). Despite 
the early comprehension of the basic emission mechanism, their spectral energy distribution is generally described by phenomeno- 
logical or simplified models. 

Aims. We propose a study of the spectral emission from the high mass X-ray binary pulsar 4U 0115+63 by means of thermal and 
bulk Comptonization models based on the physical properties of such objects. 

Methods. For this purpose, we analyze the BeppoSAX data in the energy range 0.7-100 keV of the 1999 giant outburst, 12 days after 
the maximum. We focus first on the phase averaged emission, and then on the phase resolved spectra, by modeling the system using 
a two-component continuum. 

Results. At higher energy, above ~ 7 keV, the emission is due to thermal and bulk Comptonization of the seed photons produced by 
cyclotron cooling of the accretion column, and at lower energy, the emission is due to thermal Comptonization of a blackbody source 
in a diffuse halo close to the stellar surface. Phase resolved analysis establishes that most of the emission in the main peak comes from 
the column, while the low energy component gives a nearly constant contribution throughout the phase. 

Conclusions. From the best fit parameters, we argue that the cyclotron emission is produced ~ 1.7 km above the stellar surface, and 
escapes from the column near its base, where the absorption features are generated by the interaction with the magnetic field in the 
halo. We find that in 4U 0115+63, the observed spectrum is dominated by reprocessed cyclotron radiation, whereas in other bright 
sources with stronger magnetic fields such as Her X-l, the spectrum is dominated by reprocessed bremsstrahlung. 

Key words. X-rays: binaries, pulsars: individual: 4U 0115+63 

1. Introduction The spectral emissio n from XRBPs st i ll pre sents many puz- 

zling aspects. Recently, iBecker & Wolffl d2007l) made a major 
X-ray binary pulsars (XRBPs) were discovered more than 35 step forward by proposing a model (from here on the BW model) 
years ago with the pioneering o bservations of th e bright sources w hi c h represents a complete calculation of the X-ray spectrum 
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the magnetic field orientation (see e.g. lAra va & Harding 1999; 
lArava-Gochez & Harding||2000l;ISchonherr et alj|2007k 

In this paper we use the BW model to study an observation 
of the bright HMXB 4U 01 15+634. Hard X-ray radiation from 
this source was discovered in 1969 (Whitl ock et al.| [l989). dur- 
ing one of its giant type II outbursts (16 have been observed up 
to now), when the X-ray luminosity rises to a few 10 37 erg/s, 
two orders of magnitude above its quiescent l evel. X-ray pul- 
sation s (Ps = 3.6 s) were soon discovered (ICominskv et alj 
1 19781) . and the determination of the orbital parameter (P or b — 
24.3 d, e = 0.34, sin/ = 140.1 I t - s) was made possib le for 
the first time using the SAS data dRappaport et al.lH978l) . The 
high energy (> lOkeV) continuum of 4U 0115+63 has pre- 
viously been modeled with t he standard power-l aw modified 
by an exponential cut-off (e.g. Cobur rTet al.ll2.002h . The source 
has also been widely studied in the optical and IR bands, al- 
lowing the identifi cation of its optical counterpart, V 635 Cas 
dJohns et al.ll 1 978b, and also the de termination of its distance 7— 
8kpc dNegueruela & Okazakll2001l) . 

In the spectrum of 4U 0115+63, cyclotron harmonics 
have been observed up the fifth order dW heaton et al. 1979; 
IWhite etalJ[l98H iHeindl et al.l fl999t ISantangelo et al.ll 19991) . 
while in the other known XRBP s, only the sec ond harmonic, 
and in on e case the third dTsygankov et al.ll2006l). is sometimes 
detec ted dCoburn et al l 120021: lOrlandinil 12004 ICaballero eTail 
120071) . The luminosity of 4U 01 15+63 during outburst is compa- 
rable to that of the bright sources Her X-l and Cen X-3 studied 
by BW, but the magnetic field strength, ~ 10 12 G, implied by the 
cyclotron absorption features in the spectrum of 4U 0115+63 
is considerably lower than the field strength in the other two 
sources. As a result of the analysis presented here, we find that 
the decrease in the field strength causes a fundamental shift in 
the character of the emitted X-ray spectrum, from Comptonized 
bremsstrahlung in Her X-l and Cen X-3 to Comptonized cy- 
clotron in 4U 01 15+63. 

The relation between the optical properties of the sys- 
tem and the X-ray type II outbursts have been studied by 
I Whitlock et al.l (1 19891) . wh o suggested a 3-year qua si -perio dicity 
of the outbursts . Late r, [Negueruela & Okazakj] (1200 ll) and 
iNegueruela et al.l (2001) were able to characterize the outburst 
mechanism as a result of a disc instability due to radiative warp- 
ing. 

The GINGA observations evidenced how bot h the contin- 
uum and the line energies are phase-dependent (Mihar aet alj 
120041) . A remarkable result, obtained analyzing the RXTE 
data, is the luminosit y dependence of the cyclotron features 
dNakaiima et al.|[2Q06l) . This is interpreted as reflecting the pas- 
sage from a radiation-dominated flow with a shock above the 
stellar surface at high luminosity, to a matter-dominated flow at 
low luminosity, in which the accreting material is brought to rest 
at the stellar surface by the pressure of the gas. Additional con- 
straints on the height of the line forming r egion in 4U 01 15+63 
have been developed bv lTsygankov et all (12007). who also con- 
cluded that the height of the hard-photon emission region in the 
column decreases with decreasing luminosity. 



2. Summary of the column emission model 



Af ter providing a brief summary of the iBecker &~W olff 
d2007l) column emission model, in Sect. [3] we describe the ob- 
servations and the analysis methods. In Sect. [4] we report our re- 
sults, which are discussed in detail in Sect. [3] Finally, in Sect. |6l 
we draw our conclusions. 



Beck er & Wolffl ((2007) proposed a model for the produc- 
tion of the accretion column emission spectrum based on the 
bulk and thermal Comptonization of seed photons emitted via 
bremsstrahlung, cyclotron, and blackbody processes occurring 
in the accreting plasma. The partial differential equation govern- 
ing the Green's function yields an analytical solution that can 
be convolved with the seed photon distribution to produce an 
emergent spectrum displaying a power-law continuum with an 
exponential cutoff at high energies. They assume a cylindrical 
geometry for the column, along with a constant electron tem- 
perature and a constant magnetic field. The constant properties 
of the column are a reasonable approximation, since most of 
the relevant interactions are thought to take place in a relatively 
thin layer close to its base. The density and velocity profiles of 
the infalling matter are analytical approximations of numerically 
computed functions, and the angle-dependent cross section is 
approximated using two terms to represent the mode-averaged 
electron scattering cross sections for photons propagating either 
parallel or perpendicular to the magnetic field direction, denoted 
by (cry) and (<x ± ), respectively. 

The cyclotron and bremsstrahlung source distributions each 
vary as the square of the plasma density, but their energy de- 
pendences are quite different: the cyclotron emission is due to 
the collisional excitation and radiative decay of the first Landau 
level, and it therefore displays a monoenergetic energy depen- 
dence; the bremsstrahlung emission is due to the braking of 
electrons interacting with protons, and it therefore possesses a 
continuous energy dependence. The combination of these terms 
provides a good approx imation of the fre e-free process in very 
strong magnetic fields (Riffe rtet al.ll 19991) . It should be empha- 
sized that when the cyclotron energy is comparable to the plasma 
temperature, the cyclotron emi ssion becomes the preferred cool- 
ing channel dArons e t al. 1987). Finally, blackbody emission is 
localized at the base of the column due to the presence of the 
thermal mound; it is the only term for which an integral must be 
numerically computed to derive the emergent spectrum. 

For given values of the stellar mass M* and the stellar radius 
R*, the BW model has six fundamental free parameters, namely 
the accretion rate M, the column radius rrj, the electron temper- 
ature T e , the magnetic field strength B, the photon diffusion pa- 
rameter and the Comptonization parameter 6, with the final 
two parameters are defined by 



nronipC 



M(o- ]} o- ± ) 



1/2 



S 
4 



Jbulk 
^thermal 



(1) 



where m p denotes the proton mass, c is the speed of light, and 
ytmik and ^thermal represent the y-parameter^ for the bulk and 
thermal Comptonization processes. 

The accretion rate is constrained by the observed luminos- 
ity, under the assumption of isotropic accretion and unitary ef- 
ficiency, which might allow us to remove this parameter from 
the set. A unique solution for the remaining free parameters can 
be found by comparing the model spectrum with the observa- 
tions for a particular source. The analytical nature of the model 
makes it amenable to incorporation into standard analysis pack- 
ages such as Xspec, as we did, which facilitates observational 
testing. 



1 The y-parameters describes the average fractional energy change 
experienced by a photon before it escapes through the column walls. 
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3. Observations and data analysis. 

This work is based on the observations performed by BeppoSAX 
(see [Boella et al.lll997al for a description of the mission) on 
1999 March, 26 at 17:31:32.5 (OP6714) and ended on 1999 
March, 27 at 17:34:05.5, during the decay of a giant outburst, 
more precisely 12 days after the peak, when the flux in the 2- 
lOkeV energy range was a factor of two lo wer than the maxi- 
mum value (see Fig. 1 of iHeindl et al.lll999l) . Among the many 
publicly available observations of 4U 01 15+63 during its X-ray 
outbursts, our choice is driven by the wide energy band, good 
spectral resolution, and reliable calibration of the BeppoSAX in- 
struments, and also by the relatively long exposure of this data- 
set. Another observa tion of the same o utburst, made seven days 
earlier, is reported in Santan gelo et al.l (11999b . 

We analyze the data obtained b y the high energy na rrow field 



instruments LECS (0.7-4.0 keV. iParmar et al.l 1 1 9971). MECS 



(1.5-10.5keV. lBoellaetal.1 Il997bh. HPGSPC (7-44 keV, 



iManzo et al.lll997h and PDS(1 5-1 00 keV. iFronteraet alii 19971) 
after the processing through the standard pipeline (saxdas 
v.2.1). Background subtraction is performed using the Earth oc- 
cultation for HPGSPC, the off-source pointing for PDS, and 
the standard calibration files for MECS and LECS. The LECS 
and MECS source extraction radii are set to 8' for an opti- 
mal signal to noise ratio. The instrument exposure times are 
5043 s (LECS), 53 660.4 s (MECS), 42 495.3 s (HPGSPC), and 
48 328.9 s (PDS), with the differences due to the variations in 
the operating procedures. 



3.1. Spectral analysis 

The spectral analysis is perf ormed with the Xspec package ver- 
sion 12.3.1 dArnaudfl996l) . using our implementation of the 
analytical model for the verticall y-integrated accretion co lumn 
emission spectrum developed by Becker & Wolffl (J2007)- The 
fundamental and harmonic cyclotron resonant absorption fea- 
tures are modeled using multiplicative Gaussian functions: 



G(E) = 1 



r„ V27T 



(2) 



where N n , E cn , and <x„ denote the strength, energy, and standard 
deviation of the absorption function for the nth harmonic. The 
line width is fixed whenever it cannot be determined by the fit; 
this is always the case for the fourth and fifth harmonics, set to 
the convenient value of 4 ke V. 

The phase averaged spectrum is studied at first, then we di- 
vide the pulse period into 15 bins of equal width and analyze 
each spectrum separately to obtain the dependency of the param- 
eters on the phase. Due to the short exposure, we use the LECS 
data only in the phase averaged spectrum. All of the errors are 
reported at 90% confidence level (c.l.), i.e. Ay^ jti = 2.7 for one 
parameter of interest, unless stated differently. 

3.2. Timing analysis 

We transform the photon arrival times into the Solar sys- 
tem barycenter reference frame using the standard tool for 
BeppoSAX (baryconv), t hen to the binary sys t em reference 
frame using the ephemeris (Rapp aport et al.lll978l : lTamura et al.l 
119921) . Using the HPGSPC data, we find the pulse period 
P = 3.614246 + 0.000002 s on 1999-03-26 at 17:28:45.643 UT 
(MJD 5 1 2 63 .72830605) adop ting the phase shift method de- 
scribed by Ferrigno et al] d2007l) . I.e. we divide the observation 
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Figure 1. Pulse profiles of 4U 0115+63 in different energy 
ranges. The data are from the MECS (1. 0-8.5 keV), HPGSPC 
(8.5-28 keV) and PDS (28-95 keV) instruments on board of 
BeppoSAX. 



into 174 equally spaced time bins; in each bin, we determine the 
phase of the pulse maximum with a typical uncertainty of 5% 
and then linearly fit these values as functions of time by adapting 
iteratively the initial guess on the period. The quoted uncertainty 
is purely statistical. If no orbital correction is applied, we find a 
significant variation of the pulse period within the observation: 
P = (2.6904 + 0.0005) x 10~ 9 s/s. 

We fold the light curves in different energy ranges with 100 
phase bins (Fig. [TJ. To account for the timing systematics of 
the three instruments the pulse profiles are arbitrarily shifted to 
obtain an optimal alignment in the common energy ranges (8- 
1 1 keV for MECS and HPGSPC, 20-25 keV for HPGSPC and 
PDS): we apply a forward phase shift of 4 bins (~0.14s) to the 
MECS profiles and of 5 bins (~0.18 s) to the PDS ones keeping 
the HPGSPC profiles as reference. The reason for this instru- 
mental timing problem is currently unknown. 



4. Results 

4.1. Pulse profiles 

The shape of the pulse profiles, shown in Fig. [T] is highly vari- 
able with energy: above ~ lOkeV we can clearly distinguish a 
main peak at pulse phase ~ 0.4, and, with more attention, a sec- 
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Figure 2. Pulsed fraction as function of energy, from 1.2 to 
8 keV the data are from MECS, from 8 to 30 keV from HPGSPC, 
and from 30 to 120keV from PDS. 



Table 1. Instrument 
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Figure 3. Upper panel: the unfolded spectrum of 4U 0115+63 
described by the phenomenological model: the dashed line is 
the high-energy cut-off power law, the dot-dashed line the broad 
Gaussian, and the solid line the sum. The spectra of the four 
instruments have been scaled according to the inter-calibration 
constants of Table Q] Lower panel: residuals from the best fit 
model. 



Instrument 


Constant 


LECS 


0.814+ 0.008 


MECS 


0.970+ 0.003 


PDS 


0.857+ 0.004 



ondary shallower one at pulse phase ~ 0.7 interpreted by some 
authors as the parti ally obscured emission of the column of the 
opposite polar cap (Tsvgank ov et al.l | 20 07). Below this energy, 
the pulse shape is more complex, suggesting that the lower en- 
ergy emission is generated by a different mechanism. 

We compute the pulsed fraction, i.e. the ratio ^ m ' lx ~^ min , as a 
function of energy, where the count rates at the maximum (f max ) 
and minimum (F m i n ) of the pulses are the weighted average of 
three neighboring phase bins. We perform an energy binning, 
which ensures that the signal is at least 140 times larger than the 
noise for each pulse profile. 

As shown in Fig.|2j the pulsed fraction below 5 keV remains 
always between 20% and 30%, then it begins to raise, showing at 
~ 22 keV a sharp localized decrease due to the resonant scatter- 
ing at the second Landau level. The low S/N at high energy pre- 
vents the detection of the higher harmonics, and even the funda- 
mental cannot be clearly identified despite the high count level. 
This may be due to competing effects such as photon spawning 
and cyclotron emission. We verified that in the overlapping en- 
ergy ranges, the data of the three instruments are consistent with 
each other. 

4.2. Phase averaged spectrum 

Due to the well-known problems in the absolute calibration, 
we used inter-calibration constants between the different instru- 
ments: they are obtained from a best fit procedure setting the 
HPGSPC as reference. We verified that the results reported in 
Table [TJ are model independent (within the quoted errors) and 
consistent with the instrument team recommendations Q 



http : //heasarc . gsf c . nasa . gov/ docs/sax/abc/saxabc/saxabc . html 



4.2.1. Phenomenological model 

The aim of this work is the introduction of a physical model to 
describe the spectrum. However, we study the traditional phe- 
nomenological models for comparison with the results present 
in the literature. The fit with a simple high-energy cut-off power 
law plus absorption featu res is unsatisfac t ory (y 2 d = 2.4); we 
then adopt the model that iKlochkov et all (|2008 /applied to the 
high mass X-ray binary EXO 2030+375: a broad Gaussian emis- 
sion line plus a high-energy cut-off power-law. The absorption 
lines are described by eq. (O, photoelectric absorption by the 
model phabs in Xspec. 

The unfolded spectrum is reported in Fig. [3] and the best fit 
parameters in Table [2] the continuum is altered only locally by 
cyclotron resonances, the fundamental line is shallower than the 
second harmonic as expected from the pulsed fraction of Fig. [2] 
and the photon index is almost equal to one , as typical for sat- 
urated Comptonization dCoburn et al.l l2002). The photoelectric 
absorption is parametrized by Nn = (1.38 + 0.03) x 10 22 cm -2 , 
which is higher than the value for Galactic absorption in the 
source direction, A^H.Gai = 0-9 x 1 22 cm" 2 , (average o f the val- 
ues of iDickev & Lockmanlll990tlKalberla et aT1l2005l) . thus the 
X-ray emission seems to be attenuated also locally around the 
source. 

We derive a flux in the 1-lOOkeV band of (1.69 + 
0.02) erg s _1 cirT 2 corresponding to a luminosity Lx = 6.3 x 
10 37 ergs~' assuming a distance of 7kpc. In this model, the 
power-law component is responsible for the high- and low- 
energy portions of the spectrum, which are characterized by 
completely different pulse profiles. The "bump" below the cy- 
clotron energy is described by the broad Gaussian line, which 
does not correspond clearly to any physical process. 

4.2.2. Comptonization model 

If we apply the model proposed bv lBecker & Wolffl d2007l) . us- 
ing as an input parameter the magnetic field derived using the 
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cyclotron absorption features, then the Comptonized cyclotron 
emission has a flux about two orders of magnitude larger than 
the bremsstrahlung and blackbody components, and the result- 
ing spectrum does not reproduce the lower energy portion of 
the 4U 0115+63 spectrum. Driven by the observation that the 
pulse profiles below ~ 5keV are very different from the oth- 
ers, we therefore try to fit the spectrum below this energy with 
another component. The attempts to use a blackbody, a cut- 
off power law, or a simple power-law, all attenuated by photo- 
electric absorption, are unsuccessful, giving Xied > 2.5. Among 
the reasonably simple models, only the thermal Comptonization 
of low energy (0. 1-0.5 keV) blackbody emission, located in 
the center of a spherical cloud of hot electrons, gives a sat- 
isfactory fit to the data, as verified by the applicat ion of the 
comp ST dSunvaev & Titarchukl 1 19801) and compTT dTitarchukl 
1994) models in Xspec. We will therefore attempt to fit the broad 
band spectrum of 4U 0115+63 using a combination of the BW 
bulk+thermal Comptonization model, which gives good agree- 
ment at high energies, and a thermal Comptonization component 
at lower energy. As discussed below, we interpret the low energy 
component as emission arising in a halo of gas surrounding the 
base of the accretion column. 

From the statistical point of view, it is not possible to 
discriminate between the thermal Comptonization models for 
the low energy component. We then compute the effective ra- 
dius of the blackbody which produce s the s eed photons fol- 
lowing the method of lin 't Zand et alj (1 19991) . In the compST 
model, the temperature of the seed photons is fixed to 0. 1 keV: 
for best fit parameters of the hot electron temperature T e = 
(2.91 + 0.04) keV, the optical depth of the spherical cloud t = 
22.2 ± 0.4, and the 2-100keV flux for this spectral component of 
5.13 X 10 _9 ergcm~ 2 s _I , the blackbody radius is 550 km, which 
is clearly too high for such a system. We therefore utilize the 
compTT model rather than the compST model to analyze the low- 
energy component. 

Using the compTT model, it is possible to determine the seed 
photon temperature, To = (0.51 + 0.01)keV, if we fix the ab- 
sorption column to the Galactic value. In this scenario, the low 
energy emission is produced by an optically thick plasma halo 
(t = 24.4 ± 0.4) at temperature T c = (2.90 ± 0.04) keV, while 
for comparison the BW electron temperature in the column is 
T e = (7.94 + 0.11)keV. With these parameters (reported in 
Table |2]i, corresponding to a flux in the 0.5-100 keV band of 
5.19 x 10 -9 erg cm _2 s _1 , we obtain a blackbody effective radius 
of ~16km, which is of the right order of magnitude. Assuming 
a standard NS, the gravitational redshift at the surface is 0.3. 
Correcting for this effect, we find best-fit parameters correspond- 
ing to a blackbody with ~9 km radius. However, in the rest of the 
paper we will neglect the relativistic effects for consistency with 
the BW model. 

If we fix the magnetic field strength of the BW component 
using the value derived from the absorption lines, then we find 
there is no combination of the other parameters which can ac- 
count for the emission between ~ 7 and ~ lOkeV. We there- 
fore leave the magnetic field in the BW component as a free pa- 
rameter, which effectively separates the spatial region where the 
high-energy BW component is produced from the region where 
the cyclotron absorption features are imprinted on the emer- 
gent spectrum. The implications of this physical picture are fur- 
ther discussed below. The mass accretion M is also determined 
through the spectral fit. Since there is a very strong correlation 
between the column radius ro and M, we study the contour plot 
of ro and M and find that a good fit is obtained for ro - 600 m 
and M 0.6 x 10 17 g/s. The value obtained for the accretion rate 
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Figure 4. Upper panel: the unfolded spectrum of 4U 0115+63 
described by the Comptonization model: the dashed line is the 
thermal Comptonization, the dot-dashed line the column emis- 
sion computed using the BW model, the dot-dot-dashed line a 
narrow Gaussian line, and the solid line the sum. The spectra 
of the four instruments have been scaled according to the inter- 
calibration constants of TableQ] Lower panel: residuals from the 
best fit model. 



using the fit is somewhat lower than that expected if the emission 
were isotropic, in which case the X-ray luminosity would imply 
M = 3.3 x 10 17 g/s. Our result for M therefore suggests that the 
emission from the column is anisotropic, as expected for a fan 
beam. This point is further discussed in Sect. [5] 

We then freeze ro and M to obtain a stable fit for the other 
model parameters, and the neutron star mass and radius are fixed 
at the standard values of 1.4 M Q and 10 km, which cannot be 
constrained by the comparison to the data. Among the results 
reported in Table [2] we note that the magnetic field which pro- 
duces the cyclotron emission in the BW model has a strength 
B = (5.91 +0.02) x 10 11 G while the line centroid energy implies 
that the magnetic field in the line forming region is ~ 10 I2 G. 
These results clearly suggest the existence of a spatial offset be- 
tween the BW emission region and the cyclotron absorption re- 
gion, as further explored in Sect. [5] 

Due to some residuals in the transition region between the 
low energy thermal Comptonization component and the B W col- 
umn emission, we introduce a third emission component in the 
form of a Gaussian line with energy ~ 9 keV and standard devia- 
tion <x = 1 keV. From Fig.|4]it is clear that this component gives 
a negligible contribution to the flux, unlike in the phenomeno- 
logical model, and its necessity here probably stems from our 
rough description of the cyclotron scattering in strong magnetic 
fields. 

4.3. Phase resolved spectroscopy 

The vertically-integrated accretion column emission spectrum 
computed using the BW model is independent of the observa- 
tion angle of the observer, who is assumed to be at rest in the 
frame of the neutron star. The angle independence stems from 
the fact that the model does not incorporate the full angle and 
energy dependence of the cyclotron scattering cross section, but 
instead uses a two-component formalism based on the cross sec- 
tions (tr[|) and (crO for photons propagating either parallel or 
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Table 2. Best fit parameters of the phase averaged spectral model, if uncertainties are omitted, the parameter is frozen. 



Phenomenologic al Physical 



N H [cm-'] 


1.38+0.03 


N H [cm-'] 


0.9 






To [keV] 


0.51+0.01 






T c [keV] 

T 


2 90 +0 03 
24.4+0.4 






norm* 


0.257+0.003 


E c [keV] 
E f [keV] 


14.5+0.3 
13.5+0.2 




4 +u 2 
0.59+0.02 


r 


0.998+0.007 


B[10 I2 G] 


0.591+0.002 


K 


0.247+0.003 


T f [keV] 
r [m] 

M[10 17 gs-'] 


7.94+0.11 

600 

0.6 


Eg [keV] 


7 056" 014 

' •"- MJ _n Q29 


£ [keV] 


9.07+0.12 


(T G [keV] 


2 51 + ™ 


0"g [keV] 


1 


No [ctscirr 2 s -1 ] 


0.144+0.003 


iVo [ctscirT 2 s" 1 ] 


0.0080+0.0009 


Ei [keV] 


11.16+0.15 


Ei [keV] 


11.58+0.13 


a- 1 [keV] 
N\ [ctskeV] 


1 

0.09+0.02 


O"! [keV] 
ATi [ctskeV] 


2 

0.70+0.04 


E 2 [keV] 


22.99+0.05 


E 2 [keV] 


22.76+0.09 


0- 2 [keV] 


2.8+0.1 


0-2 [keV] 


2 


N 2 [cts keV] 


2.3+0.1 


N 2 [ctskeV] 


1.22+0.05 


E 3 [keV] 


32.6+0.4 


E 3 [keV] 


36.3+0.6 


<r 3 [keV] 
N 3 [cts keV] 


3 

1.6+0.2 


cr 3 [keV] 
N 3 [ctskeV] 


3 

2.4+0.2 


E 4 [keV] 


40.8+0.7 


£4 [keV] 


46.5+1.1 


o- 4 [keV] 
N 4 [cts keV] 


4 

2.5+0.2 


cr 3 [keV] 
N 4 [ctskeV] 


4 

2.9+0.4 


E 5 [keV] 


53+1 


Ei [keV] 


57+2 


0-5 [keV] 


4 


0-3 [keV] 


4 


N 5 [cts keV] 


2.6+0.4 


N s [ctskeV] 


3.6+0.7 


t (d.o.f.) 
flux" 


1.48(705) 
(1.070+0.002) x 10~ 8 


X 1 (d.o.f.) 
flux" 


1.49(703) 
(1.070+0.003) x 10~ 8 



* the definition of the model normalization is reported in Xspec documentation. 
** the flux in the 1-100 keV energy band is expressed in units of ergcrrr 2 s -1 . 



perpendicular to the field. Despite the simplified treatment of the 
scattering cross section, it is interesting to see what insight can 
be gained by fitting the BW model to the phase-resolved spec- 
tral data. Physically, variations in the phase-resolved spectrum 
arise due to the partial occultation of the base of the accretion 
column, in combination with the angle dependence of the cy- 
clotron cross section. These physical effects can be studied in 
an approximate manner by allowing the parameters in the BW 
model to vary as functions of the pulse phase. For example, the 
effect of the changing angle of observation due to the pulsar spin 
can be modeled by allowing the cross sections (cry) and (<x ± ) to 
vary, resulting in variation of the parameter £ as a function of 
the pulse phase (see eq. (TJ). Conversely, the effect of the par- 
tial occultation of the base of the column can be modeled by 
allowing for a phase dependency in the parameters B, T e , and 5. 
Moreover, in each phase bin we note that light bending (see e.g. 
IKraus et"af. 2003) might cause the mixing of emission arriving 
from different angles. Based on these observations, we extract 
phase-dependent parameters by applying the BW model to the 
spectra in 15 equally spaced phase bins. 

For the BW component, we adopt the values of M and ro 
derived for the phase averaged spectrum. For the lower energy 
compTT component, since we do not use the LECS data, the tem- 
perature of the seed photons cannot be constrained; we therefore 
freeze it to the value 0.5 1 keV found in the phase averaged spec- 
trum. Depending on the particular characteristics of each single 
spectrum, some parameters cannot be constrained and thus are 



frozen at a suitable value consistent with the neighboring bins or 
the phase averaged spectrum. In particular, we put great care in 
limiting the line width by freezing the cr parameter in the fitting 
process whenever its unconstrained value would give a major 
contribution to the continuum model. 

We also allow a fine tuning of the relative normalization be- 
tween the different instruments due to the uncertainty in the ab- 
solute timing for BeppoSAX (see Sect. 0. However, we veri- 
fied that the scatter is within 8% (PDS) and 3% (MECS) from 
the values computed for the phase averaged spectrum and re- 
ported in Table Q] Eventually, we obtained for the fits reduced 
X 1 between 0.94 and 1.2 for 346-352 degrees of freedom cor- 
responding to a null hypothesis probability always lower than 
0.5%. Given the model uncertainties, we consider this an accept- 
able result. The unfolded spectra are reported in Fig. [5] 

The phase dependency of the best fit parameters in the BW 
component is shown in Fig. [6] the magnetic field B is higher at 
the edges of the main peak (phase 0. 13-0.66), as is the contribu- 
tion of bulk with respect to thermal Comptonization described 
by the parameter 5. The parameter ^ follows roughly the lumi- 
nosity phase dependency, whilst the temperature of the electrons 
T e displays a single peak at phase ~ 0.5. The pattern is less clear 
outside the main peak and will not be discussed further. 

The phase dependency of the optical depth t and electron 
temperature Tc in the lower energy compTT component is shown 
in Fig. |7] the former is always very high, between 20 and 35, 
and the latter varies almost by a factor of three. Even though the 
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10 100 
Energy [keV] 



Figure 5. Unfolded phase resolved spectra. The phase range 
is indicated in each plot. The dashed line is the thermal 
Comptonization, the dot-dashed line the column emission com- 
puted using the BW model, the dot-dot-dashed line a narrow 
Gaussian line, and the solid line the sum. 




0.5 1.0 

Phose 



Figure 6. Best fit parameters for the Becker and Wolff spectral 
component as function of the phase. The accretion column radius 
is fixed at 600 m and the mass accretion rate at 6 x 10 16 g s _1 . The 
solid line is the smoothed pulse profile in the 10-20keV energy 
band with arbitrary normalization. 




Phase 

Figure 7. Best fit parameters for the soft Comptonization spec- 
tral component (compTT) as function of the phase. The temper- 
ature of the seed photons is fixed to 0.5 1 ke V. The solid line is 
the smoothed pulse profile in the 10-20keV energy band with 
arbitrary normalization. 



deviation from a constant value is highly significant for both pa- 
rameters, there are no clear correlations with the pulse structure. 
The additive Gaussian line has always a centroid energy around 
9 keV with no obvious correlation to the magnetic field intensity 
(Fig. [8]), however its minimal presence is required to obtain an 
acceptable^ 2 in six phase bins. 

We show in Fig.|9]the centroid energies, the width of the first 
three harmonics, the line normalization, and the ratio between 
the centroid energies and half the centroid energy of the sec- 
ond harmonic. This final parameter can be considered the best 
indicator of the magnetic field intensity in the line forming re- 
gion, since the second harmonic is purely due to scattering and, 
incidentally, it is covered with high S/N by two instruments in 
this observation. We note that along the main peak (phase 0.12- 
0.66), the higher harmonics are almost equally spaced: the de- 
viations are at most of 4% with less than 3<x significance. Note 
that for the same phase range, the fundamental is systematically 
higher, and from phase 0.66 to 1.12, the lines do not present a 
regular spacing. 



4.4. Search for higher harmonics 

At phase 0.33-0.46 (the peak), we introduce a previously un- 
detected sixth harmonic to model some residuals in the high- 
energy tail of the spectrum. To have a better grasp on the line 
significance and verify it is not due to an inappropriate model- 
ing of the continuum, we test its presence using a cut-off power 
model for the energy range abo ve 1 3 keV. This model was used 
for previous analyses (see e.g. Santa ngelo et al.|[l999l) . The re- 
sults are reported in Table [3] the probability that the feature is 
only a statistical fluctuation is not negligible, therefore we can- 
not claim a firm detection of the sixth harmonic, just a strong 
suggestion of its presence since the centroid energy is compati- 
ble with the value of the second harmonic, the line normalization 
is greater than zero at 90% c.l., and it is present in two distinct 
phase bins at the pulse maximum. 
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Table 3. Fit results in the search for the 6 th harmonic using the cut-off power law model for the continuum above 13keV. Ay 2 is 
referred to the introduction of the further Gaussian absorption line. E a \, s is the centroid value, Eg is the value corresponding to three 
times the centroid energy of the second harmonic, and A^abs is the line normalization. The null hypothesis probability is obtained 
using the F-test. 



phase 


0.33-0.40 


0.40-0.46 


£ abs [keV] 


71+3 


68.6 ± 2.6 


E B [keV] 


66.9 ± 0.3 


70.74 + 0.15 


N abs [keV] 


4.5 ± 2.5 


3.1 + 1-3 


A^d.o.f.) 


221(168) -> 215(166) 


158(142) -> 152(140) 


Null hypothesis probability 


10% 


5% 
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Figure 9. Absorption line parameters, (a) The centroid energy, (b) The width, (c) The normalization, (d) The ratio between the 
centroid energy and half the centroid energy of the second harmonic. In each plot the 10-20keV pulse profile is sketched with 
arbitrary normalization for illustrative purpose. Line width and normalization are vertically displaced of lOnkeV, where n is the 
harmonic order. Uncertainties are reported at la c.l. assuming a Gaussian symmetric distribution. 



5. Discussion 

We have applied the self- consistent bulk+the rmal 
Comptonization model developed by iBecker & Wolffl (I2007I) to 
study the bright transient source 4U 0115+63. In order to treat 
this source, we find that the BW model needs some adjustments 



with respect to the prescriptions in the original paper. The 
first issue is the relatively low accretion rate compared to the 
X-ray luminosity: assuming a distance of 7kpc, we find a 
1-100 keV luminosity of 6.3 x 10 37 ergs _1 which corresponds 
to M = 3.3 x 10 17 gs~' for unitary conversion efficiency 
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20 



Figure 8. Best fit parameters for the Gaussian emission line as 
function of the phase. The line width is fixed to 1 keV. The solid 
line is the smoothed pulse profile in the 10-20keV energy band 
with arbitrary normalization. 



from gravitational to electromagnetic energy. If we consider 
only the BW component, we find Lx = 5.5 x 10 37 ergs _1 and 
M = 1.7 x 



10 17 gs -1 . Our adopted value for M is instead 
0.6 x 10 17 g s _1 , a factor three lower. However, using the relation 
Lx = GM„M IR„, one assumes an isotropic emission over An 
steradians. This clearly overestimates the luminosity from a 
single magnetic pole by at least a factor of two. The factor could 
be even larger if the emission came from a fan beam geometry 
and/or the NS has a different mass to radius ratio. In general, 
the prescription for the calculation of the luminosity in X-ray 
pulsars probably needs to be re-examined and firmed up due to 
these expected anisotropies in the emission profile. 

The most striking observational feature of 4U 0115+63 is 
the unique presence of at least five cyclotron lines in its spec- 
trum; other HMXBs displa y at most a second, and in one 



case a third harmonic (see Coburn et al. 2002; Orlandini 2004; 



Tsyg ankov et al.ll2006l; ICaballero et al.ll2007l) . This is due to the 
relatively low intensity of its magnetic field, the lowest among 
the Galactic HMXBs that display cyclotron absorption features. 
In fact, at energies corresponding to the higher harmonics, there 
is still a high photon density, up-scattered in energy by the ther- 
mal population of electrons at a temperature of a few keV. The 
presence of the significant high-energy continuum allows the 
harmonic absorption features to be detected. 

This peculiar characteristic finds some correspondence also 
in the continuum model. For example, we show for the first 
time that the high energy (> lOkeV) emission of 4U 0115+63 
is dominated by Co mptonized cyclotr on cooling, as expected 
on theoretical basis (lArons et al.lll987l) . Conversely, in Her X-l 
and other X-ray pulsars with higher magnetic field, the emis- 
sion is dominated by the C omptonized bremsstrahlung com- 
ponent dBecker & Wolfil 120071) . Using the best fit parameters 
obtained from the phase averaged spectrum, we simulated the 
source emission if the magnetic field were higher and found that 
the relative contribution of the cyclotron component decreases 
relative to bremsstrahlung and fades to a negligible value for 
B > 5 x 10 12 G. The reason is the decreased production rate for 
the cyclotron seed photons that form as a result of radiative de- 
excitation from the first excited Landau level, since this is less 
and less populated as its energy becomes higher than the elec- 
tron temperature. As a consequence, the spectrum of this source 
is dominated by reprocessed cyclotron emission, in contrast to 
other bright pulsars with stronger magnetic fields, in which the 
emission is dominated by reprocessed bremsstrahlung emission 
dBecker & Wolffll2007h . 

Decomposing the contributions of the different components 
as functions of the phase demonstrates that along the peak, be- 
tween 65% and 90% of the radiation is produced by the BW 



o 
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0.4 0.6 

Phase 

Figure 10. Contribution of the BW model plotted in red (dark 
gray), compTT in green (light gray), and the Gaussian emission 
line in black to the total flux (points with error bars). 



component, while at other phases the relative contribution of the 
thermal Comptonization reaches ~50% (Fig.fTOjl. The contribu- 
tion of the Gaussian emission line is always less than 5% and 
is concentrated on the descending part of the main peak. This 
is probably due to incomplete modeling of the cyclotron scat- 
tering process. The high energy component clearly originates in 
the accretion column while the low energy component has the 
characteristics of a diffuse halo, since its absolute flux is roughly 
constant throughout the phase. 

5.1. Parameter constraints 

During the determination of the optimal mass accretion rate, we 
noted that this parameter is tightly anti-correlated with the col- 
umn radius, and both are linked to £ via eq. (Q}. Since there 
is no basis for choosing a priori the value of it was neces- 
sary to study the probability contour plots of ro and M in the 
range 450m < r < 1100m and 0.4 < M u < 1.1, where 
Mn is the mass accretion rate in units of 10 17 g/s. We found 
that at 68% c.L, any combination of these parameters is ac- 
ceptable provided it falls between the lines of the equations 
r () = 1250 (M 17 - 0.56) + 450 and r = 1823 (M 17 - 0.46) + 450, 
where ro is expressed in meters. We then made the choice of fix- 
ing Mn = 0.6 and ro = 600 m, the latter according to the indica- 
tio n of a relatively l arge column in other bright sources studied 
by Becker & Wolff (2007), and also based on the upper limit se t 
bv lLamb et al] d 1973b . and later revised bv lHarding et all d!984l) . 

We also verified that the NS mass and radius cannot be con- 
strained by the spectral data. Any change in their value within 
reasonable boundaries can be compensated by tuning the mass 
accretion rate and column radius. Therefore, only four out of the 
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namely six free parameters can be determined by a direct fit of 
the model to the data. On the contrary, we already showed that 
we can constrain just a relatively wide region of the r^-M plane, 
and this region depends on the particular choice of the NS base 
parameters. Among the free parameters, the magnetic field is in- 
sensitive to the choice of the fixed parameters, while 6 and 
T e depend on them. At least in the case of 4U 01 15+63, where 
the lower energy portion of the spectrum is dominated by an- 
other component, it is not possible to infer from the data all the 
physical characteristics of the accretion column. Nevertheless, 
the BW model gives us the unique chance to investigate the ba- 
sic emission mechanisms. 



5.2. Magnetic field structure and emission geometry 

The other main issue is the strength of the magnetic field pro- 
ducing the cyclotron emission, which is found to be (0.60 ± 
0.06) x 10 12 G, while in the line forming region the field strength 
is 0.989 x 10 12 G, as derived from the centroid energy of the sec- 
ond harmonic (the values are given neglecting the gravitational 
redshift, which acts as a scale factor). Assuming a dipolar field, 
and supposing that the cyclotron absorption features form close 
to the neutron star surface, we derive the cyclotron emission is 
concentrated at ~ 1.7 km above the stellar surface. Based on this 
observation, we argue that the cyclotron cooling in 4U 01 15+63 
takes place between the NS and the column sonic point, which is 
at ~2.7km above the surface, according to the best fit model pa- 
rameters. At that height, the emission is strongly beamed since 
the plasma speed is v/c ~ 0.25 and most of the radiation is emit- 
ted downwards, as predicted for an emittin g charge in relativis- 
tic m otion (see e.g. eq. 4.100 - 4.103 of iRvbicki & LightmarJ 
fl979h . 

The radiation is subject to Compton reprocessing inside the 
accretion column, and it can subsequently escape from the sides 
only near the base, where the advection velocity decreases, and 
the density and temperature are higher. We can verify that for 
the best fit value of the diffusive escape time from the col- 
umn is compar able to the free-fall tim e from the sonic point - 
see eq. (104) of lBecker & Wo lff (2007). Further confirmation of 
this scenario comes from the phase dependency of the parame- 
ters £ and T e , which are higher along the peak, implying respec- 
tively that the radiation traveled mostly parallel to the magnetic 
field, and that the visible layer is deeper in the column for this 
phase range. The sharp feature in the pulsed fraction (Fig. |2]i 
at ~ 20 ke V also fits into the picture: while the emerging radia- 
tion is beamed roughly perpendicular to the column, the resonant 
scattering produces an isotropic emission which reduces locally 
the pulsed fraction. 

The lower energy component is due to thermal 
Comptonization of seed photons with temperature 
Tq = (0.51 ± 0.01)keV, which is not far from the temper- 
ature of the thermal mound probably present at the base of the 
column: according to the BW model of the high energy emission 
we find Tjh ~ 0.7 keV. However, any attempt to fit the broad 
band spectrum using the BW model alone was unsuccessful 
due to the very low flux of the components originating from 
the blackbody and bremsstrahlung sources compared to the 
cyclotron one. As was already discussed, the effective radius 
of the blackbody photon source for the low energy compTT 
component is about 15 km (or 9 km considering the gravitational 
redshift). It is therefore likely that the low-energy radiation 
source resides close to the NS surface, perhaps in a diffuse halo. 
An optically thick atmosphere (t ~ 20) at a temperature of a 
few keV could be confined by multipolar or crustal components 
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Figure 11. A schematic and not in scale cartoon of the magneto- 
spheric structure of 4U 01 15+63, as it emerges from our analy- 
sis. 



of the magnetic field (see e.g. iGil et al.ll2002l) . while the region 
around the column could be heated to ~ 0.5 keV either by the 
accretion itself or by the strongly beamed column emission. 

The modulation of the magnetic field producing the contin- 
uum radiation along the main peak (Fig. |6]l suggests that the 
emission layer has not a planar geometry, but becomes closer 
to the NS at the edges of the column. A similar pattern is present 
in the phase dependency of the 6 parameter, indicating that bulk 
Comptonization is stronger at the center of the main peak. The 
photon-diffusion parameter, follows nicely the pulse profile in 
the phase interval 0.2-0.8. This can give us some information on 
the angle between the magnetic field and the radiation direction. 

Even though £ is inversely proportional to the geometrical 
mean of the extreme values of the electron scattering cross sec- 
tion (see eq. IT]), and does not include explicitly its angle depen- 
dence, it can be regarded as an estimator of the average value in a 
particular phase bin. We argue that at the maximum of the peak, 
the radiation field is subject to a reduced scattering cross sec- 
tion, and therefore the photons must t ravel preferent ially in the 
direction of the magnetic field lines dVenturalll979l) . The pho- 
tons are therefore able to reach almost the NS surface, before 
escaping from the column walls mostly near the base, where the 
plasma temperature is higher. At the column boundaries, where 
the plasma is optically thin outside the cyclotron resonances, 
the observed cyclotron absorption features are imprinted on the 
spectral energy distribution. In Fig.QT]we provide a schematic 
drawing of the model that we propose. 

5.3. Standard spectral phenomenological model 

Even though we were able to successfully model the continuum 
with a sta ndard spectral phenom enological model similar to the 
one that iKlochkov et ail J20081) used for EXO 2030+375, we 
consider this model unsatisfactory for 4U 0115+63 due to its 
incoherency with the pulse energy dependency. The pulse pro- 
files and the pulsed fraction change at ~ 7 keV, corresponding 
to the spectral region where the Gaussian "bump" starts to dom- 
inate, whilst they do not exhibit a peculiar evolution at the en- 
ergy where the power law continuum becomes again the main 
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component. It is also difficult to explain why the same spectral 
component produces a pulsed fraction which is nearly constant 
below ~ 5 keV and steadily increasing above ~ 15 keV. 

Adopting the spectral phenomenological model, we note the 
appearance of some residuals in the spectral fit below ~ 1 keV 
(Fig. |3). Even though they are not statistically significant due 
to the low exposure of the LECS instrument, they essentially 
disappear when we adopt the multi-component Comptonization 
model (Fig.|4j». With this latter model, we find that the photoelec- 
tric absorption due to neutral plasma is entirely due to the fore- 
ground Galactic medium: this is plausible since the observation 
took place after the outburst maximum when the strong radiation 
field of the NS could have ionized the surrounding medium. 

Moreover, to our knowledge there is no physical process 
which could explain the broad Gaussian "bump" in emission: 
for instance, if it is due to cyclotron cooling, then one must find 
a way to make the Compton scattering so ineffective in repro- 
cessing this component that it maintains a broad Gaussian pro- 
file instead of the typical power-law with exponential decay. We 
conclude that although this model can yield a good phenomeno- 
logical description of the spectrum, it is nonetheless unable to 
give a coherent picture of the system in terms of physical pro- 
cesses. 

The decoupling between the high and low energy emission is 
more naturally explained by two distinct emission components, 
one dominating below ~ 7 keV, and the other above. In the sce- 
nario we have focused on here, the low energy emission (48% of 
the phase averaged total flux in the 1-100 keV energy band) is 
modeled as a th ermal Comptonization of 0.51 keV seed photons 
dTitar chuk 1994), while the high energy continu um (5 1 % of the 
total fl ux) is modeled using the BW formalism (Becker & Wolff 
120071) . We also add a Gaussian emission line which gives a mi- 
nor contribution to the total flux (~ 1%). The appearance of this 
component likely reflects either an approximation of the emis- 
sion wings predicted to appear near the fundamental by Monte 
Carlo simulations of cyclotron scattering (Schon herr et al.l l2007. 
and references therein), and/or an adjustment to the BW model 
to account for a more complex geometry in the scattering pro- 
cess. While this component is required to obtain an acceptable 
X 2 , it is clearly not a major contributor to the continuum, in con- 
trast to the phenomenological model described above, in which 
it contributes 16% of the total flux in the 1-100 keV energy band. 



5.4. Structure of the cyclotron harmonics 

The line widths of the second and third harmonics reach a min- 
imum just after the maximum of the pulse, which indicates that 
the line of sight reaches its maximum angle relative to the mag- 
netic field at this phase. However, we note that the lines ap- 
pear much shallower t han those computed in the MC models 
(Schonher r et al.ll2007l) . Probable reasons are a more complex 
geometry with respect to the model assumptions and general rel- 
ativistic effects such as light bending, lensing effect, and doppler 
shift. A more quantitative discussion of this issue and the sce- 
nario we envisage for the system would require a complete nu- 
merical simulation, but this goes beyond the scope of the present 
paper. 

If the radiation field dominates the thermodynamics of the 
plasma, then one expects a link between the temperature of the 
plasma and the magnetic field intensity. The relation depends 
on the optical depth in the absorption lines, the angle that the 
radiation form s with the magnetic field lines, and the system ge- 
ometry (e.g. Ilsenberg et al. 11998, and references therein). In a 
slab geometry and for radiation injected parallel to the magnetic 



field, the cyclotron energy should be more than two times the 
plasma temperature. For isotropic radiation, one instead expects 
a ratio of about 4. This is the case for the halo temperature and 
the cyclotron line energy found here, while the column temper- 
ature is too high to be due to cyclotron heating, supporting the 
idea that the line forming region is not located within the column 
but rather in an external halo. 

The apparent non-harmonicity of the line spacing in some 
HMXBs has been explained b y the variation of the magnetic 
field in the line forming region. Nishimura (2003) showed that 
for a dipolar field the height of the region relevant for the ab- 
sorption and scattering process must be of the order of several 
hundred meters t o reproduce the observations. In a later work 
iNishi mura (2005) demonstrates that a linearly varying magnetic 
field is able to produce line non-harmonicity even for a 5 m thick 
slab of plasma above the NS. Such a configuration can be mod- 
eled by a combination of axial and crustal magnetic field, as pro- 
posed by [GUetaD (Hooi). 

We show that the higher har monics are equal l y spac ed in 
4U 0115+63, as also noted by ISantangelo et al.1 d 19991) and 
Tsyganko v et al.l (|2007), and therefore the field is nearly con- 
stant in the line forming region. The previously reported devia- 
tions from harmonicity regards only the ratio to the fundamen- 
tal, whose centroid energy can be displaced from the cyclotron 
value. This can be due to the influence of the other emitting com- 
ponent, which is present below the cyclotron energy, and/or to a 
complex shape of the line profile, heavily modified by photon 
spawning, and then poorly described by a Gaussian absorption 
line. Also the need to introduce a Gaussian emission line points 
towards a complex structure of the cyclotron fundamental line, 
with emission wings that are not modeled by a Gaussian profile. 

Using the second harmonic, with energy ~ 20keV, as the 
best estimator of the magnetic field in the line forming region, 
the displacement from the harmonic law is always less than 5% 
for the higher harmonics, thus the scattering region thickness is 
at most 2% of its distance from the NS center. The modulation of 
the line centroid ene rgy with phase was already reported e.g. by 
Miha ra et al.l ([2004 ) and is a typical feature of this source: it is 
interpreted as due to the occultation of part of the column by the 
NS along its rotation. Since the magnetic field of the line forming 
region varies from ( 1 .07 + 0.08) X 1 12 G to (0.87 ± 0. 1 5) X 1 12 G 
along the main peak, the maximum reduction from the maximum 
value is 23%. This corresponds to an 8% variation of the height 
of the line forming region; assuming that at phase 0.5-0.6 the 
absorption takes place close to the base of the column and the 
NS has a radius of 10 km it means that for instance at phase 0.2- 
0.3 the scattering takes place in a layer, which is a few hundred 
meters thick, located ~ 800 m above the NS surface. 

The last issue we would like to briefly discuss is the puzzling 
jump of the energy of the fundamental cyclotron feature around 
the pulse minimum from ~ 11 to ~ 15keV. A similar discon- 
tinuity was a l ready n oted in the pha s e aver a ged spectrum by 
Miha ra et al.l (120041) ; iNakaiima et all d2006l) : iTsygankov et al.l 
(120071) when the luminosity drops below 5 x 10 37 erg/s: it is in- 
terpreted as a drift of the compact emitting region in the mag- 
netic field, due either to the disruption of the radiative shock 
or possibly to the accretion column burning out with the on- 
set of another column in a multi-polar configuration. We show 
that this jump is not linked exclusively to the system lumi- 
nosity, but appears at different phases even when the source 
is brighter than the proposed threshold for the shock disrup- 
tion (~ 6 x 10 37 ergs _1 ). Moreover, the absorption features at 
30^4-0 keV which are present outside the main peak, do not 
seem to be linked by an harmonic ratio to the ~ 20keV line, 
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which suggests that complex field geometries could be present 
in 4U 01 15+63. 



6. Conclusions 

For the first time we apply the self-consistent model of the 
high luminosity X-ray binary pulsar emission proposed by 
iBecker & Wolff! (120071) to the well studied transient source 
4U 0115+63 during the giant outburst that occurred in 1999. 
We find that the emission of the peak is produced by thermal 
and bulk Comptonization of the cyclotron emission which is the 
main cooling channel of the column. On the other hand the lower 
energy continuum is roughly constant throughout the phase, and 
is due to thermal Comptonization of a 0.5 keV blackbody, whose 
effective radius is comparable to the NS dimensions. 

We are able to sketch a model in which the photons are gen- 
erated by cyclotron cooling well above the NS surface. The pho- 
tons are advected towards the NS by the relativistic plasma bulk 
motion, while they are Compton up-scattered by the hot elec- 
trons in the column. The radiation escapes from the column lat- 
eral walls close to the NS surface, where the relativistic beaming 
is less severe. At this height the cyclotron lines are formed, prob- 
ably in the boundary layer between the column and the plasma 
which forms an extended halo above the NS surface, and also 
produces the low-energy spectral component. 

Unequivocal interpretation of the phase-dependent parame- 
ters using the theoretical model would require some estimate of 
the angles for the system. I.e., the angle between the line of sight 
and the rotation axis, and the angle between the rotation axis 
and the accretion column. These two angles would determine 
the portion of the column visible during a given phase, as well 
as the angle between the normal to the column and the line of site 
to the observer. Relativistic ray tracing should be used to study 
the influence of the system geometry on the emission pattern, 
but this is beyond the scope of this paper. 

Although we have developed a convincing physical model 
for 4U 01 15+63 which is able to describe coherently the system 
properties, we emphasize that this scenario should be checked 
quantitatively in future work using a detailed simulation. The 
analysis presented here was mad e possible via the imp lemen- 
tation of the model proposed by Becker & Wolfl] (120071) in the 
standard fitting package Xspec. The application of this new 
model is very promising and it is expected to facilitate the study 
of many other systems. 
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